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Within the Dirac model for the electronic excitations of graphene, we calculate the full polarization 
tensor with finite mass and chemical potential. It has, besides the (OO)-component, a second form 
factor, which must be accounted for. We obtain explicit formulas for both form factors and for the 
reflection coefficients. Using these, we discuss the regions in the momentum-frequency plane where 
plasmons may exist and give numeric solutions for the plasmon dispersion relations. It turns out 
that plasmons exist for both, TE and TM polarizations over the whole range of the ratio of mass 
to chemical potential, except for zero chemical potential, where only a TE plasmon exists. 


I. INTRODUCTION 

At present, graphene is still an object of highly actual interest, especially for its optical properties and plasmonics. 
There is an enormous number of papers on, for recent reviews see [1] and [2]. One of the approaches to the theoretical 
description of its electronic excitations rests on the Random Phase Approximation (RPA) for the density-density 
correlation function and for the conductivity. Usually, these quantities are calculated in some approximation, the 
non-dispersive limit for example. Another approximation rests on the smallness of the Fermi velocity, uf, which is 
300 times smaller than the speed of light. As shown below, this smallness allows formally to neglect one of the two 
form factors entering the polarization tensor. 

The electronic properties of graphene are equally well described by the Dirac model, consisting of a relativistic 
spinor in a (2+l)-dimensional space-time (the plane of graphene). It has its own Lorentz group with up in place of 
c. The coupling to the (3+l)-dimensional electromagnetic field is the usual one. In this way, one has a modified 
Quantum Electrodynamics (QED), to which the well-known formalisms can be applied. For instance, one calculates 
the photon polarization tensor IRA In the language of Quantum Field Theory (QFT) it consists of a fermion loop 
(if restricting to one loop approximation). This has been done in quite a number of papers, for instance in [3] and 
related papers (mainly in application to external magnetic field). For a kind of review we refer to [4]. In fact, this 
approach is in principle equivalent to RPA, where, however, one makes some approximations for the very beginning. 
The formulas obtained using the polarization tensor allow to calculate the Casimir force between graphene and, say a 
conducting wall, ideal [5] or real [6], and to investigate surface plasmons, e.g. in the TE mode [7]. Surface plasmons 
were considered earlier using expressions for the conductivity obtained in RPA. For instance, in [8], using a high 
frequency approximation for the conductivity, a TE mode plasmon was predicted. 

In the present paper, the approach from QFT is taken. The complete polarization tensor is calculated for both, 
mass m (gap parameter) and chemical potential p, non-zero. In this case, like with finite temperature, the polarization 
tensor, which is transverse for gauge invariance, has two form factors, and both are accounted for. Also, the complete 
frequency and momentum dependence is kept. All approximations, usually done within the RPA approach can be 
obtained afterwards as special cases. 

The polarization tensor was previously calculated in a number of special cases. First of all one needs to mention 
that the calculations in the RPA approach deliver special cases for components of the polarization tensor. So in [9, 10], 
in fact the (OO)-component of the polarization tensor, which is responsible for static screening and for the TE mode 
plasmon, was calculated in the massless case. In the QFT approach, this component was calculated with chemical 
potential and mass in [11]. 

In application to graphene, the second form factor was calculated in [12] with mass, temperature and chemical 
potential, however restricted to Matsubara values of its argument, i.e., to discrete imaginary frequencies. It was applied 
to the calculation of the Casimir effect. Recently, it was calculated for finite temperature and mass for all, including 
non-Matsubara and real, frequencies [13], but without chemical potential. From these papers, a representation with 
mass and chemical potential for real frequencies, as it is needed for plasmons, cannot be derived in any simple way. 
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Here we fill this gap and calculate the complete polarization tensor for both non-zero, mass and chemical potential, 
for real frequencies. 

We apply this polarization tensor to the investigation of surface plasmons. For this we start from the generic 
formulas for the electromagnetic field in (3+l)-dinrensions and insert the polarization tensor resulting from the (2+1)- 
dimensional fermions. We show that this polarization tensor, which involves the Fermi speed uf, can be expressed 
in terms of the same polarization tensor calculated with unit speed. The actual calculation is done in the Appendix. 
Further we (re-)derive the equations for the reflection coefficients in terms of the two form factors, Eq. (32), and, 
which is equivalent, in terms of noo and n tr , Eq. (38). Further, we solve the equations for the plasmon dispersions 
numerically and represent the results graphically for the TE and TM-nrodes. 

Throughout the paper we use units with a = e 2 /( 47 r) for the coupling and v = Vp/c for the Fernri speed and refer 
to 


1 1 

137’ V ~ 300’ 


(1) 


as the physical values of the parameters. 


II. ELECTRODYNAMICS WITH POLARIZATION TENSOR FROM GRAPHENE 

By the Dirac equation model, the interaction of the electromagnetic field with the long wavelength electronic 
excitations in graphene can be described by relativistic quantum electrodynamics. In this formalism, the Dirac 
equation is 


where 


(i$ — eA — m)ijj = 0 , 


r?P = diag(l, i>, v), 


with the usual gamma matrices 7 ^ and the Fermi speed v. The resulting polarization tensor is 


=ie 2 


dq 0 d 2 q 1 




1 


(27t) 3 ij — m itj[ — if — m 


( 2 ) 

( 3 ) 

( 4 ) 


where we used momentum representation. Assuming the graphene sheet located at z = 0, the momenta are in the 
directions ju = 0,1,2. It can be reduced to the polarization tensor \V iv (p ), defined by the same formula but with 
v = 1 , by the substitution qi -+ qi/v (i = 1 , 2 ) in the integration, 

n^) = ^<' nMV (M', ( 5 ) 

where the vector p is substituted by p. This vector is defined by 

■ ( 6 ) 

We will use the notation with a tilde for vectors with the Fermi speed v in the spatial components and for the 
polarization tensor (4) throughout the paper. 

In the four dimensional formulation of electrodynamics, the Maxwell equations with polarization tensor (which we 
denote by a ’hat’), read 


d^ v + fl^A^ = 0. (7) 

These can be viewed as effective Maxwell equations, which appear, e.g., integrating out the spinor fields in a functional 
integral representation. Rewritten in terms of induced current, 

f = ^a v , ( 8 ) 

47T 

these equations are = (4 n/c)^. 
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Switching to 3-dimensional notations with /i = (0, k) (A: = 1,2,3), 

f = ( cp,j ), A* = (M), E = ^V4> - d 0 A, 

we get from the transversality of the polarization tensor 

= fl^d^Ek, 

and the induced charge density and current, 

p = n“V^. 

A = 

We would like to mention that from this formula the conductivity tensor can be defined as (Jki 
The Maxwell equations (5) read now 

div£? = 47rp, 

A'jr 

{-d 2 0 +A-VoV)£= — d 0 j, 

C 

or, substituting Gauss’s law, 

( -d^ + A) E = (\7cp + d 0 fj . 

Inserting from (11), the right hand side can be expressed in terms of the polarization tensor, 

(-dl + a) E k = - (v fe n ol do 1 + n fci ) e,. 

These are general formulas. For graphene we have to insert 

frv = ( 6(z)iU“'(p) for = 0,1,2, 

| 0 for p, = 3 or v = 3, 


(9) 

( 10 ) 


( 11 ) 

(-c/A-n^do 1 . 


( 12 ) 


(13) 


(14) 


(15) 


with W v {p) from (4). 

In momentum representation we assume all relevant quantities ~ exp (—iuit + *kx) and define the vectors in the 
plane of graphene by 



and similar for other vectors. Further we need to define the vector 

p^ = (uj, k, 0), p = \/u} 2 — k 2 . (17) 

Now we split equations (13) into two with n,m = 1, 2, i.e., parallel to the plane of graphene, and and a third in the 
perpendicular direction, 

(.P 2 + d 2 z ) E n (p, z ) = S(z)E nm E m (p, 0), 

(p 2 + d 2 z ) E 3 (p,z) = d z S(z)-U On (p)E n ( P ,0 ), (18) 

10J 

where we introduced the notation 

= —n 0m (p) - n nm (p), (19) 

LO 

appearing in the right side in (14). In (18), we have to solve the equations for the components E n (p , z) of the electric 
field parallel to the plane, whereas E 3 (p,z) follows from E n (p, z) by integration. 
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Now we need a more specific expression for the polarization tensor for graphene. From (5) we get 

n 00 (p) = 4 n °°(p)> 

ri tr = ^—^n 00 (p) + n tr (p), (20) 

where we defined 

n trip) = n(p) tr = g^(p) (21) 


for the traces. 

The polarization tensor is calculated in the Appendix for v = 1. There it is represented in terms of form factors 
A(p) and B(p), Eqs. (A9), where we have to insert with ( 6 ) 


= (w, vk, 0 ), p= \/u 2 — (uk ) 2 


( 22 ) 


in the form factors and in the tensor structures, (A10). For these we note 


V 2 k n kn 


P Um (p) = -frS P nm (p ) = s*—, 

M 0m (p) = M nm (p) = 


Accounting also for the we get from (5) 


—wkn 


p z 


jjOm _ 

fl nm = - [ S r , 

This allows to rewrite (19) in the form 


m o~B(p) 




v 2 k n kr, 

p 2 


+ fWhl B(p). 

pi 


= s nm - (1 - ir) 


2 , k n km 


A(p) - 


»V k n k v 

9 ~9 

L0 Z P Z 


B{p). 


(23) 


(24) 


(25) 


The form factors A{p) and B(p) are calculated in the Appendix, Eqs. (A12) and (A53), (A65). 

Next we introduce the polarizations for the electric field. We consider only the components E n with n = 1,2, 

E n (p, z) = ^ 'Fte + ^tm, (26) 

where, up to a normalization, we have Ate = $te and Atm = 9 2 $tm- We do not need to consider here the third 
component, A 3 , since it follows with (18) from A„. Now, Eq. (18) is diagonal in the polarizations, i.e., the graphene 
does not mix these. So we get with (26) from (18) 

(. V 2 + )$t x(p, z) = d(2:)5 T x < l ) Tx(p, 0). (27) 


Here the subscript 'TX' stands for one of the polarizations and 


1 (-hi 


-TE = 


-TM = 


k 2 V h 
1 (h 


M n 


-fo 

h 


= A(p), 


k 2 \k 


M n 


h 


= A (p) -v 2 ~2 B (p) 


UT 


Finally, we rewrite the equation (27) for z = 0 in terms of matching conditions, 

d>Tx(z = +0) — d> T x(z = —0) = 0, 

d z $Tx(z = +0) — d z $Tx(z = +0) = S T x < 1 ) tx(2 = 0), 


(28) 


( 29 ) 



5 


for both polarizations, i.e., the functions are continuous and their derivatives jump. In a standard scattering setup, 
the solutions are 


$tx(*) = (e ipz + r T xe~ ipz ) Q(-z) + t T xe ipz <d(z) 
and the reflection and transmission coefficients are 


^TX 


-1 

1 + Q t -^ 


Itx 


1 

1 + Qtx ’ 


with 

QtE = 2^ A(P) ' Q ™ = W ( MP) l ’ 2 ' 

These are the final formulas for the coefficients. Similar formulas were derived earlier, e.g., Eq. (23) in [12]. 


(30) 


(31) 


(32) 


III. SURFACE PLASMONS 
A. General formulas 

Surface plasnrons appear if the reflection and transmission coefficients (31) have a pole. So these are solutions of 
the equations 


1 + Q tx = 0. 


We look for them in the frequency region 


vk < lo < k. 


(33) 

(34) 


The upper bound implies a frequency below the continuous spectrum in order to get a wave function (30) decreasing 
to both sides of the graphene sheet. In this region the momentum p is imaginary and we use the notation 

p = ir] = %\J—u) 2 + k 2 (35) 


and get from (32) 

Qte = ^A(p), Qtm = ^ ^A(p) - v 2 —B(p)^j . 

Here we insert from (A12), 

A (p> = + n “(p)> 

These form factors can be inserted into (36) and we get 

Qte = ^ (ilp-n 00 (p) + n tr (p)) = i(gn 00 (p) + n tr (p)), 
Qtm = jvhc 7 n 00 (p) = 2^n 00 (p), 


(36) 


(37) 


(38) 


where we also displayed the expressions in terms of the IP 1 "^), Eq. (4), using (20). Eqs. (38) are the final formulas 
for the investigation of the equations (33). The components of the polarization tensor n M "(p) are given by Eqs. (A54) 
and (A65). These coincide with Eq.(24) in [12], or Eq. (12) in [13], where, however, different notations are used. 

It is meaningful to check the case p < m, which we considered in [7]. In that case we have B — 0 for the second 
form factor and from (A19) the relation 


— 2 p 2 

v 2 k 2 


n 00 (p) 


n tr (p) 


(39) 
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holds. From (37) we get 


m = 


In the notations of [7] the polarization tensor is expressed in terms of the function 


by 


delivering 


2 / fj 

$(p) = - 2 mp — (p 2 + 4m 2 )arctanh-f— 
p \ 2 TO 


n ™<« = -SF IjF ■"<« 


= Tr, 5 <5™ = ‘5 5 


This coincides with Eq. (18) in [7]. 


(40) 


(41) 


(42) 


(43) 


B. Formulas for the surface plasmon in the massless case 


Frequently, graphene is considered with zero gap widths, which seems justified given the small value of the gap 
widths. In the Dirac model this translates into the massless case. For m = 0, the formulas (31) and (38) remain 
unchanged whereas those for the polarization tensor simplify. We collected them in the Appendix B. With the 
substitution p->pwe get from (B2) and (38) 


Q te — 
Qtm = 


4a 
77 

4 ap 
(vk) 2 


u > 2 p pi 
(vk) 2 


4 2 


sign(< 2 ) ^ x\/x 2 — 1 + arccosh(a;)^ 


Ai=±l 


(vk) 2 1 

4 p 2 


sign(Q) ^ x\/x 2 — 1 — arccosh(a;)^ 


Ai=±l 


The variables are 


77 


= \Jk 2 — w 2 , p = \/u> 2 — (vk) 2 , x = ———^—, Q = p 2 + 2 Ai up. 


vk 


(44) 


(45) 


As said above, surface plasmons appear as solutions of the equation (33). Now, inserting numbers using (50) shows 
that for the TE polarizations the solutions start with k > 0, whereas for the TM case the solutions start from k = 0. 
In this case, the solution can be considered for small k. A direct expansion of (44) in powers of k gives 


Qtm = —a 


Vk 2 — uj 
u> 


2 , 2 u + w . 

In -z— -b 0(k). 


2/7 — U) 


(46) 


This expression is in agreement with Eq. ( 8 ) in [ 8 ], which turns out to be the approximation for small k and uj. Eq. 
(33) with (46) can be solved for small k by iteration. First we rewrite the equation in the form 


u> 


= o.\Jk 2 — uj 2 In — — + 0(k). 

2/7 — UJ 


Next we insert uj = 0 in the right side and get cj = cj s f(k) in the left side with 

x>sf(k) = ak + 0(k 2 ). 


(47) 


(48) 


This solution confirms that we get small ui s f(k) for small k as anticipated in the derivation of (53). It is interesting to 
note, that this solution, because of vk < w s f (k) < k, (34), is consistent for 


v < a < 1 

only, which is fulfilled for the physical values of the parameters. 


( 49 ) 
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C. Regions in the (fe,w)-plane 

As said in Eq. (34), plasmons may exist in a region in the (fc,w)-plane below the continuous spectrum. In addition, 
the polarization tensor must be real. Since we are in the Minkowskian region, it is real below the threshold at 

to s = \j ( vk) 2 + (2 ?n) 2 , (50) 

i.e., for u) < u> s . In the case without chemical potential, i.e., for p = 0, this is the only region where the polarization 
tensor is real. With chemical potential there is also above the threshold a region of realness. It is bounded from above 
by the condition (A68), which reads 


UJ < U) + 


(51) 


with 


oQ = fi + y/p 2 + (vk) 2 — 2 vk kF (52) 

written in the notations used here. Combining (34) and (51), we get 

vk < ui < min(fc,u;+), (53) 

which defines the region in the (k, w)-plane, where plasmons may exist. 

This region is shown in Fig. I 1 . It looks different for different ratios to/ p. For to = 0 it is shown in Fig. 1(a). Here 
degenerates, = p + \p — vk |, and the region (53) becomes a triangle. Since here io s = vk holds, all solutions 
are automatically above the threshold. 

For in > 0 we introduce the following notations, 


k s 


2m 

\/l — v 2 ’ 


2 [n — vkp) 
1-v 2 ’ 


k t 


2 kp 
v ’ 


(54) 


where fci? = \/n 2 — m 2 is the Fermi momentum. In (54), k s is the intersection k = u> s , k m is the intersection k = u + , 
and kt is the touching point = u) s . For to = rn t , with 


m t 


\/1 — V 2 [l, 


(55) 


these coincide, 


k s = k m = k t = 2/j (to = TO t ). (56) 

The case to < nit is shown in Fig. 1 (b) and mt < to < /x is shown in Fig. 1 (d). It must be mentioned that the last 
region is very narrow, 


M - m t 


5-10- 6 , 


and TO t = 0.999994/r for physical values of the parameters. 


(57) 


D. TM mode surface plasmons 

TM mode surface plasmons are solutions of the equation (33), 

1 + Qtm = 0, (58) 

where Qtm is given by (A54) and n 00 (p) by (A58) with p inserted for p. 

We start with the case to = 0. Here, Qtm is given by (50). As said above, the solution exists in the triangular 
region shown in Fig. 1(a). Examples are shown in Fig. 2, left panel, for several values of a and v = 1/300. These 


1 All plots are made for p = 1. 
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FIG. 1. The (fc, cu)-plane with the curves for vk, k, cu 3 and ca7 for (a) m = 0, (b) m = 0.6, (c) m = 0.999984/r, i.e., m < mt, 
and (d) m = 0.999997, i.e., mt < m < /r, for physical values of the parameters. The region (53), where a plasmon may exist, 
is shaded. In panel (b), kt. km and outside the graph. The curves u> p and u> s appear as straight lines and the line vk nearly 
coincides with the fc-axis. The inset shows the same picture for larger k such that kt is seen. Now the line c o = k nearly 
coincides with the to axis. Panel (d) corresponds to the narrow region mt < m < /i, where the touching point kt. is below k s 
and km- 



FIG. 2. Left: Transverse magnetic plasmon with m = 0, fi ^ 0, and v = 1/300, for different values of a. The inset shows 
K = — lo%- m as a function of vk/^ for a = 1/137. Right: Transverse electric plasmon. These plots are very close one to 

the other differing significantly only by their endpoints. The inset (a) is the zoom of TE curves. The nearly horizontal dashed 
line is u> = 2/r — vk. The inset (b) shows \fk 2 — lo^, e as a function of vk/[i. 


solutions start in k = 0 and terminates on the line w = 2^ — vk. Since these solutions are all restricted to k < v//i, 
the spatial extend of the solutions in the sense of exp(— *c\z\) is determined by m: = \Jk 2 — ui 2 . It is shown in Fig. 2 
in the inset. For the physical value a = 1/137, the TM solution bends closer to vk. A similar picture may be found 
in [14]. 

For m > 0, the plasmon solution exists in the regions shown in Fig. 1, (b)-(d). The curve w7 separates from ui = vk 
and in the gap between them a new branch of the solution appears for large k, as shown in Fig. 3, panel (a). 

Obviously this branch and the lower one are parts of a single solution which, however in the gap between these 
branches is not real. When further increasing m, starting from m = mo, these branches merge, see Fig. 3, panel (b). 
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The mass, for which the solution touches the curve oj + at its minimum, k m i n = \J /. i 2 — m 2 / 2, is denoted by mg. It 
can be found by substituting k m i n and u)+(k m i n ) = m + n, into Eq. (58). With physical values of the parameters, 
this equation yields mo = 0.23235p. 



FIG. 3. TM-plasmon dispersion law for physical values of the parameters. Panel (a) shows the solution with m — 0.1 < mo 
consisting of two branches. Panel (b) presents the plot for m = mo = 0.232347, corresponding to the mass when two branches 
of the solution merge. For this mass fc max = 597.863. Panel (c) with m = 0.6 demonstrates that with increasing mass the 
endpoint moves toward the origin. The curves for ojT, uj s and vk (from top to bottom) are shown as dashed lines. 


The solution has an endpoint, 

; y/3^/8a 2 p(p, — m) (1 — v 2 ) + 3 m 2 v 2 3m 

max = Wi^ Wi^’ ( } 

which moves down to smaller k when further increasing m, see Fig. 3, panel (c). Finally, for m —> p, i.e. when the 
chemical potential disappears, this endpoint goes down to k = 0, and the solution disappears. 

For small k it is possible to find the solution explicitly. For this, we rewrite Eq. (58) in the form 


u= \jk 2 - 4u 2 fc 4 IIo 0 2 , 


(60) 


which allows for iteration with inserting ui = k in the right side, 

„ 4 i 


= k- 


v 


8a 2 /z 2 (1 — 


0(k 4 ). 


y/m 2 v 2 +n 2 (l-v 2 ) ' 


(61) 


The coefficient in front of k 3 is negative, showing that the solution goes indeed below ui = k, i.e., below the border 
of the continuous spectrum. This coefficient is small, proportional to (v 2 /a) 2 ~ 2.3 • 10 -6 , unless m —> ji, where it 
becomes infinite and the expansion breaks down (expansion (48) holds instead). 


E. TE mode surface plasmon 

TE mode surface plasmons are solutions of the equation (33), 

1 + Qte = 0, (62) 

where Qte is given by (38) with (A54) and (A65) with p inserted for p. The TE mode surface plasmon solution has 
a nonzero staring point, which we denote by ho- It can be found as a solution of the equation 1 + Qte|^ =! , = 0. For 
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k < ho, the solution goes into the continuous spectrum. It is shown in Fig. 4, together with k s and k rn . Eq. (54), as 
a function of to. The intersection of these curves in the left panel is denoted by Too- This means, for to > too, that 
the starting point in the (k, w)-plane, see Fig. 1, is to the left of k s and for to < too to the right. The latter implies 
that the starting point is above the threshold. 

In the case m = 0, as mentioned above, the solution exists in the triangular region shown in Fig. 1, panel (a). 
Examples are shown in Fig. 2, right panel, for several values of a and v = 1/300. The TE solution goes close to u> = k. 
Its frequency lies in the band 1.667 < w/n < 1.9934. This is in agreement with [8]. The starting point lies at oj = k 
and does not depend on a. The endpoint is situated on the line ui = 2/r — vk and depends on a , see the inset. With 
decreasing a it tends to k m = 2 /x/(l + v). 

For in < m t , the solution ends on the line oh/, thus above the threshold. Such solutions are shown in Fig. 5 for 
several m < m t . In the ( k , w)-plane, all appear nearly on one and the same line. Therefore we represented them in the 
(k, w//cp)-plane. Close to the end, these lines have a knee. This is shown in Fig. 5 in the right panel for a particular 
value of to using smaller values of a enlarging this part. 

Further increasing to, for to > mt, the endpoint appears below the threshold and it goes to larger k the closer to 
comes to /i. At the same time the starting point goes to zero, see Fig. 4. For to —> fi, which is the transition to the 
case without chemical potential, the starting point is zero and the endpoint goes to infinity. The picture turns into 
that shown on Fig. 2 in [7]. 


ko 



FIG. 4. The starting point ko of the TE-plasmon as a function of m for physical values of a and v, together with k s and k m . 


O) 


at 



2 P 
2m 




a> + 



m = 0.9 


0 12 3 4 5 
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k 


FIG. 5. Left panel: TE-plasmon dispersion laws for physical values of the parameters (solid lines) for several values of the 
mass, m = 0,0.4, 0.6, 0.8, 0.9, as function of k/kF- The uppermost curve corresponds to zero mass and is the same as in Fig. 
2, right panel. Right panel: the same for m = 0.9, for several values of the coupling, a = 1,0.5,0.2,1/137. The knee at the 
end of the curves in the left panel becomes visible only for coupling stronger than its physical value. The curves for oj^_ and 
LOs are shown as dashed lines. These are nearly horizontal. Also the straight lines for c o = k are shown as dashed lines. The 
corresponding line for vk coincides in the plot with the fc-axis. 


IV. CONCLUSIONS 

In the foregoing sections we calculated the polarization tensor for graphene with mass and chemical potential using 
formalism and notations of quantum field theory. 
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The actual calculation of the polarization tensor is carried out in the Appendix, for unit speed of light. In the case 
considered here, all integrations can be carried out and explicit formulas for the complete polarization tensor emerge. 
It is seen that it does not have an ultraviolet divergence, which is the expected result from the interplay of gauge 
invariance and dimensionality. 

The polarization tensor with to and /i has 2 form factors. The result is given in terms of these. Formulas of this 
kind for the polarization tensor were obtained previously, but never in such completeness as here. Especially, in [11], 
only one form factor, II 00 (p), was calculated. 

It should be mentioned, that the second form factor, II tr (p), Eq.(A13), frequently gives only a small contribution. 
This can be seen from the upper line in Eq. (38), where Iltr(p) enters with a relative weight proportional to v ~ 1/300 
as compared to II 00 (p). In other cases, its contribution is essential. For example, for p < in, from (39), n tr (p) gives 
twice the contribution as II 00 (p) and has opposite sign. Thus dropping n tr (p) in Qte, Eq. (38), changes even the 
sign, which results in the disappearing of the TE plasmon in this case. 

The explicit form of the final formulas looks different in different regions in the (fc, w)-plane. All these are related 
by corresponding analytic continuation, which can be best understood in terms of the left side of Eq. (A37), viewed 
as an integral representation of the result. For the regions relevant for plasmons, i.e, where the polarization tensor is 
real, we gave explicit expressions in terms of real functions, Eqs. (A54), (A65) for frequencies below the threshold, 
and Eqs. (A70), (A71) for above. Also, we gave in Appendix B explicit formulas for the massless case. 

In section II we started from the general 4-dinrensional notations of QED, used the effective Maxwell equations and 
related these to the (3+l)-dimensional formulation with (2+l)-dimensional polarization tensor. This allowed, after 
separation of the polarizations into TE and TM, to formulate a scattering setup for the electromagnetic field and to 
relate the scattering coefficients with the form factors of the polarization tensor. We represented the whole derivation 
in necessary detail to make the paper self contained. Also, at the end of the section, we discussed the restoration of 
v ^ 1, including the mixing of the from factors, Eq. (20). Finally, in Sec. III.A, from the reflection coefficients, we 
derived the equations for the plasmons, including the simplifications appearing for in = 0, Sec. III.B. 

In the remaining part of Sec. Ill, we investigated the plasmons. There, a quite sophisticated number of cases 
appeared. First of all, the regions in the ( k , w)-plane, where plasmons exist, quite strongly depend on in in the 
interval 0 < m < p, see Fig. 1. In the massless case, we reproduce known results, see Fig. 2. When the mass starts 
to grow, these pictures deform, see Figs. 3 and 5. Finally, for to = /./, i.e., when the chemical potential disappears, 
we turn into the case considered in [7]. 

A general property of these solutions is that there are always (except for p = 0) solutions for both polarizations 
present. Thereby, for small in, the TM solution is present starting from to ~ k > 0 and the TE solution only for 
larger k and w, see Fig. 2. For to closer to p, the range of k , where the TM solution exists, shrinks, Fig. 3, panel (c), 
whereas the TE solution extends to larger ^-intervals, see Fig. 5, right panel. Roughly speaking, for to = 0, /z ^ 0, 
there is a TM solution and for m ^ 0, p = 0, there is a TE solution. 

It must be mentioned that, possibly, not all solutions considered here, are physical in the sense that some, e.g., Fig. 
3, left panel, exist only for k exceeding the range of validity of the Dirac model for graphene. Another observation 
is that the smallness of the physical parameters, enters in different way. Besides the general proportionality to a of 
the polarization tensor, and to the Fermi speed v in the range (34), the combination v 2 /a enters, Eq. (61). The 
combination a/v, which is of order of unity, did not play any pronounced role. 

It would be interesting to carry on a similar analysis of the plasmons for finite temperature without and with 
chemical potential. 
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Appendix A: Polarization tensor in (2 + 1) dimensions with chemical potential 

In this appendix we display the calculation of the polarization tensor. We consider a space-time with metric 
= diag(l, — 1, — 1) with Greek indices p, v = 0.1.2. For the spatial part we use Latin indices i,j = 1,2. As said 
in Section II, a metric with a Fermi speed wp / 1 can be restored afterwards. Following [5, 12, 15], the polarization 
tensor is 


B M1/ (p) = *2e 2 J tr S(q)^S(q - p)j‘ 


(Al) 



12 


with the spinor propagator 


S(q) = 


id + to 


q 2 — to 2 — i 0 


(A2) 


and for the 4-dimensional gamma matrices { 7 ^, 7 1 '} = 2 g^ v holds. The factor 2 in front comes from the number 
Nf = 2 of flavors. In this appendix vectors are denoted by 


V M = (Po,p), P = (PUP 2 ), 


p = \/pI-p 2 , p = \Jp\ + pi¬ 


rn 


In this appendix we use lower indices for the time components, IIoo = II 00 , for convenience. The units are taken such 
that a = e 2 / 47 r holds. Carrying out the trace in (Al) we get 




(A4) 


where 


= <f(q v ~ Pi + (<f - Pld" - q(q ~ p)<T + ™ 2 g 
N = {q 2 — to 2 + iO)((q — p) 2 — m 2 + iO), (A5) 

with q M = (q 0 , q), q = |q|. For instance, we note 

^00 = qo(qo -po) + q(q - p) +to 2 . (A6 ) 

The go-integration in (Al) is specified by the causal '+ iO'-description. The polarization tensor (Al) is transversal, 

P^ V (P) = 0. (A7) 


Before carrying out the integration this manifests itself in the structure of the numerator which can be written as 

P^ = -((q- p) 2 - m 2 )q v + {q 2 - m 2 )(q - p) v . (A 8 ) 

In each term a factor from the denominator (A5) cancels and the remaining integrals correspond to tadpole diagrams 
and vanish due to parity. 

In case of chemical potential or temperature, the polarization tensor depends, besides on also on an additional 
vector nP — (1,0,0). From transversality we have 2 independent tensor structures, 


I= P^(p)A(p) + M^(p)B(p), 


(A9) 


with 


pn p v 


P^{p) = g^ M^(p) = 


pfipV p\l n v _|_ n fipV 


np 




w) 2 


Below we will need the special cases 


Poo(p) = -2*, M 00 (p) = 

g^P^ip) = 2, g^M^( P ) = -g. 


(A10) 


(All) 


The form factors can be obtained from (A9), 

A{p) = oo(p) + n tr (p), 

pZ 

B(p) = 2^n 00 (p) + 4 n t rip), (A12) 

P 4 P Z 


where 


n tr (p) = g tlv W‘ v (p) 


(A13) 
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is the trace over the Lorentz indices. We will use this subscript also for other tensors. Further we note 

n°" = ^n 00) (A14) 

p z 

which is a direct consequence of the transversality (A7). 

For the calculation of n tr we mention with (A5) 

Zti = -q(q - p) + 3m 2 

= -^(g 2 - m 2 + {q - p) 2 - rn 2 ) + P 2 + 4?n 2 ), (A15) 

where in the second line the g-dependent terms are written in a way cancelling the corresponding factors in the 
denominator, and a term, which is independent of q. In this way we get 



n t ,(p) = -s tp + p2 + 2 4m2 s(p), 

(A16) 

where 

v, . 0 2 [ dq 0 d 2 q 1 

tp J (2 tt) 3 g 2 — m 2 + *0 

(A17) 

corresponds to 

a tadpole graph and 



- lie 2 f dqod2q 1 

yP> J (2 t r) 3 (g 2 -m 2 + i0)((g-p) 2 -m 2 +*0) 

(A18) 


is the scalar loop. 

We mention that in case there is no vector n ,J in , the tensor structure is given by the first term in (A9) alone 
and the relation 


n 00 (p) 


h^ {p) 


(A19) 


must hold. In this way, the calculation of the full polarization tensor can be reduced to the calculation of n 0 o or of 
the scalar loop E and the tadpole E tp . 

Next we calculate IIoo and II tr directly. In the calculation we consider IIoo, E and E tp in parallel since most steps in 
the calculation are the same for them. Starting from here we include the chemical potential p. For this we include for 
a moment temperature too using the Matsubara representation. Thus we substitute po by the Euclidean momentum, 
po = ip 4 with p 4 = 27 tIT (1 integer) and we turn the go-integration to the imaginary axis and include the chemical 
potential, go = *g 4 — p with g 4 = 27r(n + 1/2 )T (n-integer). Indicating p in the argument we get 


n 00 (p-,n) 


-8e 2 Tj2 

n 


d 2 q (iq 4 - p)(iq4 - ip4 ~ p) + q(q - P) + m 2 
(27r) 2 [(ig 4 - p) 2 - Tf] [(*g 4 - ip4 - p) 2 - r|] 


St P (p;p) = -8e 2 T^ 
n ' 

E(p; M ) = - 8 e 2 T^ 


d 2 q 


1 


(2 tt ) 2 (m-p) 2 -^ 
d 2 q 


(27r) 2 [{iq4 - p) 2 - rj] [(ig 4 - ip4 - p) 2 - r|] 


with the notations 


r 4 = y/q 2 + m 2 , r 2 = y/(q — p) 2 + m 2 . 


(A20) 


(A21) 


We mention, that IIoo in (A20) coincides up to an overall factor with Eq. (A.3) in [11]. 

Further in this paper we restrict ourselves to F = 0. So we go back from the Matsubara summation to the 
integration, Tj2 n —t / dg 4 /(27r). To proceed, we rewrite IIoo(p;/4) and E (p; p) in a form with denominators linear in 
g4, 

<%d 2 g 1 M\ 1 \ 2 


n 0 o (p\p) = Se 2 J 


E 


(27r) 3 4 ^ [ig 4 - p + AiTi] [ig 4 - ip4 — p + A 2 T 2 ] ’ 

Ai,A2=±1 


Z(p\p) = —8e 2 


dq^d 2 q 1 ^ A 1 A 2 


1 

p + AiTi] [ig 4 - ip4 - p + A 2 T 2 ] ’ 


(A22) 
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with 


M AlA2 =l + ^i^-(q(q-p) + m 2 ). (A23) 

t it 2 

This rewriting is done in such a way, that the momentum (74 appears in the denominator only. 

Now, in order to proceed, we need to handle the ultraviolet divergence. In general, there are none in the polarization 
tensor in (2+1) dimensions. This is obvious from power counting in (A20). We have 3 powers from the integration 
and -2 powers from the propagators. Due to the gauge invariance, 2 divergent powers drop out and a convergent 
expression is left. However, in intermediate steps, before the compensation due to transversality is on work, there 
are divergences. Therefore we assume a regularization not breaking the gauge invariance and not affecting the ( 74 - 
integration. Obviously this is possible. An example is dimensional regularization in the spatial directions. 

With such regularization assumed we carry out the ^-integration in (A22) using 



dx 

(ix + a)(ix + b) 


sgn(i?e(a)) — sgn (Re(b)) 

71 -;- 

—a + b 


which follows simply with the Cauchy theorem. With a —> 


—71 + AiTi and b —> —ip 4 — /i + A 2 ^ we get 


n 0 o(p;M) = - 8 e 2 


d 2 g 1 y- s gn (—71 + AjTi) - sgn(-p + A 2 T 2 ) 

(2tt) 2 8 AlAa — AiTi — ip4 + A 2 F 2 

Ai ,A 2=±1 


E(p;/z) = - 8 e 2 


d 2 q 1 y- AiA 2 sgn(- 7 i + AiTi) -sgn(-/x + A 2 r 2 ) 

(27r) 2 8 f- TiTs -A^ - ip A + A 2 T 2 

Ai ,A 2=±1 


The signs factor in the numerators can be rewritten as 


sgn(-7i + Aifi) — sgn (—p + A 2 T 2 ) — Ai — A 2 — 2 (Ai©(Ai 7 t — Ti) — A 20 (A 2 /i — T 2 )) 


(A24) 


(A25) 


(A26) 


using the step function. Now, for p < m, the expression in the parentheses vanishes and only the contribution without 
chemical potential is left. This allows to separate the contributions from the chemical potential, 

n 0 o(p; m) = n oo } (p) + A M n oo(p), 

E(p;m) = E ( 0 ) (p) + A m E(p), (A27) 

where the superindex ’( 0 )’ denotes the vacuum part and 


A M n 00 (p) = — 8 e 2 


A m E ( p ) = — 8 e 2 


d 2 q 1 
(27r) 2 4 

d 2 q 1 
(27r) 2 4 


y. a7 Ai a 2 

Ai ,A2 =il 


AiQ(Ai/i — Ti) — A20(A2 /z — T 2 )) 

ip4 + AiTi — A2T2 


E 

Ai ,A 2 =il 


Ai A 2 

fy 


Ai0(Ai/i — Ti) — A20(A2 /x — T 2 ) 

ip4 + AiTi — A2T2 


(A28) 


are the additional contributions from the chemical potential. Similar separation takes place for the tadpole contribution 
which, however, will be calculated later directly. 

In the above expression, the momentum integration is bounded by Ti < \p\, T 2 < \p\ and the integration is finite. 
This is in line with the general situation that in the polarization tensor in an external field (e.g., a magnetic field) or 
with temperature can be separated into the vacuum part and the field or temperature dependent part, which is free 
of ultraviolet divergences (at least in one loop). 

Once the momentum integration is bounded, we can remove any regularization and calculate the integrals directly. 
For this we exchange in the second term in the numerators Ti -O- T 2 , Ai O A 2 . The exchange of the T’s appears from 
the substitution <7 —» p — q of the integration variable in (A28). We mention that this is in general not possible if a 
regularization, a momentum cut-off for example, is present. From (A28) we get this way 


A M n oo(p) = ^ 8 e 2 


d 2 q 


4 E M Xl x 2 Q(Xip-T 1 ) 

Ai ,A2=±1 


—Ai 

ip4 — AiTi + A 2 r 2 


Ai 

ip4 + AiTi — A 2 r 2 


^) = - 8 « 2 /(Cy 1 E 


A 1 A 2 


^ 4 A,.fc ± , r ' r » 


0(Ai/i — Ti) 


—Ai 


Ai 


ip4 — AiFi + A 2 r 2 7+4 + AiTi — A 2 T 2 


(A29) 
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Next, we reorder the summation by substituting A 2 = crAi, 

/ i2 1 

, 9 ? 2 4 ^2 M a Q(X iZ-i-Ti) 

X > Ai,<t=±1 L 


_ ~~Ai _|_Ai_ 

ip4 — AiTi + AictT 2 ip4 + AiTi — Aicrr 2 J 


A „E(p) = — 8 e 2 J 


2 f d 2 q 1 


< 2 ') 2 4 r > r > 


0 (Ai/x —Ti) - 


?P4 — AiTi + AioT 2 *p 4 + AiTi — AioT 2 


. (A30) 


Now the expression in the square bracket is in fact independent on Ai which takes values ±1 only. Thus the sum over 
Ai goes only over the step function and becomes 


E ©(Ai/x-r 1 ) = 0(/x 2 -r 2 ) 

Ai=il 


(again we used lb >0). We arrive at 


A M n 00 (p) = - 8 e 2 [-2-^ 0(M 2 -r 2 )i E M °~ - rr 

J (27 tY 4 ^ ip4 — AiTi 

Ai ,cr=±i 

= -** 2 / ^ e (, 2 -rl)l y ^ ^ 

Ai ,(T=±I 


+ Aio-r 2 


[_ -Ai _ 

r 2 ip4 — AiTi + Aio-r 2 


Now it is meaningful to turn to polar coordinates for the integration. For this wee need to remove the square roots 
containing the cosine which can be done using 

_El_ = r i ^?'Aip 4 + oT 2 

ip4 — AiTi + AioT 2 (ip4 — AiTi) 2 — T 2 

Carrying out the summation over a we get 


A,U 00 (p) = -8e 2 J^Q(p 2 -Tl)l E 

Ai — ±1 

a ^) = - 8e2 /(0 0( ' ,2 - r?) 5 ? SF 


z = r 1 (r 1 -iAip 4 )+q(q-p)+m 2 , 
n = ( i P 4 - A^!) 2 - r 2 . 

The angular dependence is in qp = qp cos <p and in T| = T 2 + p 2 — 2 qp cos <p and we rewrite 

Z = ^ (-JV + (2T, - *Aip 4 ) 2 - p 2 ) , 

IV = —p 2 — p 2 — 2 iAiP 4 T 1 + 2gp cos <p = Q + a cos tp. 

Now the angular integration can be carried out by the formula 

f dtp 1 1 

J 2n Q + a cos tp v /Q 2 _ a 2 ’ 


and we get 


A „n„ o( p). Ad r dri i v if ,' i ir --; a «)’-p’ 

" w ' 2jr J m 1 2 ^ 2 ^ ^ 7^2 




1 

2 _ „ 2 ’ 


Ai=±l X/Q ° 
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where dqq = oSTiTi was used. The square root is uniquely defined at least for both real, Q and a, with \Q\ > |a|. 
For complex Q as in (A36) the sign of the imaginary part of the square root follows the sign of the imaginary part of 
Q, i.e., it is — Ai. 

In fact we started the calculation of the polarization tensor in the Minkowskian region and turned to Euclidean 
momenta only temporarily. Now, after carrying out the ^-integration we turn back by substituting p 4 = —ipo- Then 
Q in (A36) becomes Q = p 2 — 2Aip 0 ri> where (A3) was used. The radicand in the denominators in (A38) can be 
rewritten, 


2 2-2 
Q -a = p p 


2Ti — XiPo\ 4m 2 — p 2 


P 


V 


with 


x = 


2Ti — Aipo 

— y 

V 


9 4m 2 

x o = — 

pZ 


_ 2 —2 / 2 I 2 \ 

= p p {x + x 0 ) 


- 1 . 


(A39) 


(A40) 


Under the rotation p 4 = —ip o, for Ai = +1, the Q moves toward the real axis. For \Q\ > a, which we assume in the 

(A41) 


following, the square root \JQ 2 — a 2 becomes real and its sign follows the sign of Q. Thus we get 

\JQ 2 - a 2 -)• sign (Q)pp\Jx 2 + Xq 

and the last square root must be taken positive. Rewriting in the first line in (A38), 


(2Ti - iAip 4 ) 2 -p 2 =p- 


we get from (A38) 


A^IIoo^) — 
A M £(p) = 


—8e 2 r1 


— AiPq 

P 


dFl i E \ (-H-sign(Q)^ 

Ai=±l \ 


- 1 


= P\x 2 - 1), 


2?r J m 

—8e 2 f M 


p ( x 2 — 1) 


P \/x 2 + Xq 


1 


27r 


^ E 


sign(Q) 


2 Ai—=ti PP\/x 2 + ■ 


Now the last integration can be carried out using the indefinite integrals 


dx 


dx 


\/ x 2 + x 
x 2 — 1 


= arcsinh- 


x 0 


\Jx 2 + x% 2 


xyj x 2 + Xq — (2 + Xo)arcsinh 


and dTi = dxp/2. 

Denoting the indefinite integrals by noo(p;Ti) and £(p;r 4 ) we get 


n 0 o(p;Ti) = E ~ x + sign(Q)^ ^ x^Jx 2 + x\ - (2 + a^arcsinh 


Ai — i 1 


X 

Xo 


£(p;Ti) = -^-^ E sign(Q) arcsinh—. 

dnp 2 x ^ ±i xo 

Taken in the corresponding boundaries these give 

A M n 00 (p) = n 00 (p; p) - n 00 (p; m), 

A M £(p) = £(p;/x) - £(p;m). 


(A42) 


(A43) 


(A44) 


(A45) 


(A46) 


We will see that with the functions (A45) just that functions appear which were defined in (A20). 
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For this to show we consider these functions for T\ = in and show that they coincide with the corresponding vacuum 
expressions in (A27). In (A45), the dependence on Ti is in the variable x, defined in (A39), only. For Ti = m , which 
is equivalent to q = 0 (see (A21)), we get x = (2m — \iPo)/P and note 

x 2 +xl = -^= 2 - (A47) 

pZpZ 

This follows with (A21) from (A36) for q = 0. We have to take the square root, 

sign(Q)^/a; 2 + x% = (A48) 

where the sign follows from (A41). Further, with T = m and pa = —ioj, we get from (A36) 

Q = p 2 - 2Ai p 0 m, (A49) 

where p 2 = Pq — p 2 from (A3). These formulas allow to express 


sign(<3)-\/ x 2 + x'q _ P 2 - 2Ai p 0 m 
x p(2m — Aipo) 


(A50) 


Further we use arcsinhz = arccoth-\/z -2 + 1 and arccotli = arctanha — arccoth)- to get 


arcsinh— = arccoth ^ + A 


xo 


= sign(Q) arccoth 


p 2 — 2Xipom 
p(2m - Xipo)' 


= arctanh^ - A^rctanh^. 

2m po 


(A51) 


Using the above relations we rewrite (A45), 
—8e 2 p 1 I" 2m — Xipo 


II oo(p;m) = 


8t r 2 


I E 


Ai==hl 


P 


1 - 


pp 2 -2X 1 p 0 m\ p 


2P 


pp 


-Z- 1 + 


) 2 p 


4 rrT 


( arctanh JL - Al arctanh^ 
2m po 


v 


E(p;m) = 


—8e 2 1 


47rp 2 


E ( arctanh-^-Aiarctanh 

V 2m 


Ai=±l x 

The summation simplifies these expressions, 


Po 


(A52) 


n 00 (p; m) = ^ P 3 (imp - {p 2 + 4m 2 ) arctanh^-) 


167Tp' 

. —8e 2 p 

E(p;m) = —-arctanh-—. 

47rp 2m 


(A53) 


Now the statement is that Eqs. (A53) give just the vacuum contributions in Eq. (A27). For both expression this can 
be verified easily, especially for the scalar contribution. For the polarization tensor we mention Eqs. (6) and (7) in 
[5] and Eq. (A. 16) in [11] (with vf = 1). 

In this way, when adding (A45) with the vacuum expression in (A27), the latter cancel and we are left with the 
final representations 


n 0 o (p; p) = 


8e 

4-7T 


|Ea 1= 


±1 

—8 e 2 p 2 
167T p : 


sign(Q) p‘ 
4 p 


[x\/x 1 + Xq — (2 + Xq) arcsinh^-^ (p > m) 


!- (2 mp — (p 2 + 4m 2 ) arctanh^-) (p < m), 


and 


Ef ) = / ^ * Ea 1= ±i sign(Q) arcsinhJL > TO ), 

I Tif arctanh^- (p<m), 


(A54) 


(A55) 
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where now 


x = 


- Aipo 


4m 2 


x n = 


— 1, Q=P + 2Xnxp. 

p p* 

It remains to calculate the tadpole contribution in (A20). We rewrite it for T = 0, 

d 2 q 1 


s t P (p;p) = Se 2 J ^ J 


(27t) 2 (iq 4 - p + ri)(*g 4 - p-T i) ’ 
where we wrote the denominator as a product. Now we apply Eq. (A24), 

d 2 q lsgn(-/z + ri)-sgn(-/z-ri) 


St P (p\p) = -*]■ 


(27r) 2 2 -2Ti 

The sign factors can be rewritten in the form 

sgn(-ju + Ti) - sgn(-p- Ti) = 2 - 20(p 2 - T 2 ), 
allowing to split into the vacuum part, 


S^(p) = -8e 2 / 


d 3 


- 2 e r 


1 


(2 tt) 3 - 2£ —g 2 - q 2 -m 2 ’ 


and the addendum 


A M £ tp (p) = 


—8e 2 /' d 2 q 0(p 2 -T 2 ) 


(A56) 


(A57) 


(A58) 


(A59) 


(A60) 


(A61) 


2 J (2t r) 2 ' r. 

In the vacuum part we introduced dimensional regularization since the integral is divergent. Direct integration gives 


_ o„2 r ( £ - l/ 2 ) _ -8e 2 
(47t) 3 / 2_£ 47T 


K’(p) = 8e 


m + 0(e). 


(A62) 


As usual, in even space dimension the dimensional regularization does not give a pole contribution. 
In (A61), there is no angular dependence in the integrand and the g-integration is trivial, 


Putting (A62) and (A63) together we get 


—8e 2 

A M S t P (p) = O - m )- 


V ~ 8e2 ( \ 

E tp = —— max(/i, to), 


47T 


(A63) 


(A64) 


which is the final result for the tadpole contribution. 

Now, from (A16) with (A62) and (A61) we get for the trace of the polarization tensor 


Htr (p-,n) = 



p+4m 1 


2 P 


l Ea 1= ±i si g n (<3) arcsinh^) , (p. > m), 


8 e 2 

8np 


(2 mp — (p 2 + 4m 2 ) arctanh^) , (p < m). 


(A65) 


By equations (A55) and (A65) we have explicit expressions for the two form factors of the polarization tensor. In the 
course of the derivation, after Eq. (A38), we turned from the Euclidean to the Minkowski region. Thereby, especially 
in Eq. (A44), we assumed Xq to be real. With (A40) this implies p 2 < (2m) 2 , i.e., the region below the threshold, 
which is, in terms of frequency at po = to s with co s = \Jp 2 + 4?n 2 . Thus, Eqs. (A55) and (A65), as written, are in 
this region in terms of real functions. 

The above formulas allow for an easy analytic continuation to the region above threshold, i.e., to po > uj s . In the 
complex po-plane, the polarization tensor has a cut starting from io s . We consider the continuation in the upper half 
plane, i.e., for Spo > 0- Thereby we have 


x 0 = 


I (2m) 2 . 

2 -=2 - 1 with 2/0 = 

Pq-P 


I Po ~ (P 2 + (2m) 2 ) 
Po - P 2 


(A66) 
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Note p 2 > p 2 because we are in the Minkowskian region. For the continuation in Eqs. (A55) and (A65) we use the 
relations 


arcsinh(iz) = = sign(;j) + arccosh(|,z|)) , 


7T / 1 

arctanh(i^) = i —I- arctanh — 
2 V 2 


(A67) 


(\z\ > 1), where we have to insert z Using (A41) we note x 2 — = ( Q 2 — a 2 )/(p 2 p 2 ), thus \x/yo\ > 1. Using 

Eqs. (A56) and (A36) (for a), Q 2 — a 2 can be seen to be a fourth order polynom in p 0 . Its roots are 


(A) = Xfj, + cr'y/p 2 + p 2 + 2ap Icf (cr, cr'= ±1) 


(A68) 


Analyzing these roots, the restriction po < w7 follows, together with p < po- Further we need the sign of Q, Eq. 
(A56). It is a second order polynom in p 0 and from analyzing its roots under the above restrictions, sign(Q) = Aj 
follows. 

Using now (A67) in the upper lines of Eqs. (A55) and (A65), and accounting for the sum over Ai, the imaginary 
contributions cancel. As a result we get 


sign(Q) arcsinh 
Ai 


—> ^ Aiarccosh 
Ai 



(A69) 


for the transition to above the threshold. In all other places in (A55) and (A65) we substitute simply Xq —> — z/q. 

In the expressions for m > /./, i.e., in the second lines of Eqs. (A55) and (A65), we use the second line in (A67). 
Together, the final formulas for the from factors of the polarization tensor above the threshold read 


and 


Hoo (p;p) 


( §£. 

J 477 

M + Ea 1= ±i A i (W^ 2 — 2/o — ( 2 — Vo) arccosh ( 

X 

yo 

))] 

1 

illfi (2 mp {p 2 + 4m 2 ) ( 1 1 + arctanh 2 ™)) 

(m < 


n tr 0; At) 


8e2 

477 

8e 2 

8 np 


(m + p ~ + p m2 h E Al =±! A 1 arccosh ( 

(2mp — (p 2 + 4m 2 ) + arctanh^pj ^ 


(p > in), 

(p < m). 


The upper lines of these equations are real for 


(A70) 


(A71) 


vV + (2m) 2 < po < uj 1 . 


(A72) 


This defines the region where the polarization tensor is real which is of relevance for the plasmons. In all other regions, 
one has either an imaginary part or, for physical values of a and v, the frequency is outside the validity region of the 
Dirac model. 

As for the lower lines, i.e., for m > p which is equivalent to no chemical potential, there is the usual imaginary part 
above threshold resulting from particle creation. 


Appendix B: The case m = 0 

Her we specify the formulas (A54) and (A65) for the case m = 0. In this case only the upper lines of these formulas 
apply. For the variable x still (A56) applies. For Xo we get x\ = —1 and taking ui with positive imaginary part, 
xq = —i. Further we need 


arcsinh(zx) 


iarcsin(a;), for \x\ < 1, 

sign(a;) + arccosh(|x|)) , for \x\ > 1 


(Bl) 


For |x| < 1 we observe imaginary parts from this formula and from the square root in (A54) as well and in this region 
we cannot expect a stable solution. In opposite, for |x| > 1, the expressions become real after the sum over Ai is 
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done, which makes the contribution from in/2 vanishing. Further, if we consider /.t > 0 only, we have x > 1 and can 
drop the signs in (Bl). In this way, in the massless case the polarization tensors become 


n 0 o (p;m) = 


2e 2 


H~ ^ sign(Q) (x\/x 2 — 1 — arccosh(a;)^ 


Ap —i 1 


2e 2 

n t r (p;ju) = — 

7T 


p ^ 

/J - 2 2 sign(Q) arccosh(x) 

Ai=±l 


(B2) 


with 


2^ - AiPq 

P 


(B3) 


which is the same as (A56). 
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